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We use analytic techniques to study the gravitational force that would be produced by different 
Kuiper-Belt mass distributions. In particular, we study the 3-dimensional rings (and wedge) whose 
densities vary as the inverse of the distance, as a constant, as the inverse-squared of the distance, 
as well as that which varies according to the Boss-Peale model. These analytic calculations yield 
physical insight into the physics of the problem. They also verify that physically viable models of 
this type can produce neither the magnitude nor the constancy of the Pioneer anomaly. 



I. INTRODUCTION 

There has long been interest in the gravitational force 
that could be produced by the Kuiper Belt [l|. It has 
been observed that total masses of much more than an 
Earth mass, M®, would lead to conflicts with orbital ob- 
servations. (See, e.g., Refs. 0,0 and Sec. VII-E of 0.) 
Further, it has also been calculated that a Kuiper-Belt 
ring with a mass of this magnitude could not explain an 
acceleration the size of the Pioneer anomaly [J Q 
This anomaly 0, E| is the apparent unmodeled constant 
acceleration of the Pioneer spacecraft, observed between 
~ 20 — 70 Astronomical Units (AU), of magnitude 

a P (20AU < r < 70AU) = -(8.74 ± 1.33) x 10~ 8 cm/s 2 

which is directed approximately towards the Sun. 

Even so, this type of Kuiper-Belt mechanism has re- 
mained a fascinating one as a possible explanation of the 
anomaly. In particular, it has recently been proposed 
[j| that gravitation from the Kuiper Belt, modeled by 
a cylindrically-symmetric ring of matter whose density 
goes as 

Pi(p) = P = V x2 +V 2 , ( 2 ) 

P 

where 

pi = 1.74 x 10- 16 g/ cm: \ L = 20 AU, (3) 
can explain the constant anomaly. The ring has a width 
Ri = 20 AU < p < 100 AU = R 2 (4) 
and a thickness 

2D = 2 AU. (5) 

The mass is thus 

Mring = 4ir Pl L D(R 2 - Ri) = 1.17 x 10 28 g = 1.96M e . 

(6) 
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This proposal is somewhat surprising, given the obser- 
vations noted above. However, one is thereby motivated 
to take a different looks at the problem . Here we do so 
emphasizing analytic calculations. This will help to bet- 
ter understand the underlying physics of the situation. 

To start, although it is well-known that a spherically 
symmetric ball with a density that goes as l/r can pro- 
duce a constant acceleration within the ball, there only 
is a constant acceleration from a complete spherical ball, 
not from a shell. Therefore, as we emphasize in the 
next section, with only a cylinder ring, not even a cylin- 
drical disk, satisfying a constant acceleration is doubly 
hard to do. Specifically, it can not come from an exact 
cylindrically-symmetric 1/p density. Indeed, although 
the appeal to Gauss' Law in Eq. (3) of 5] is correct, 
the argument that Eq. (4) of [5j implies there will be a 
constant acceleration within the ring is not exact. We 
will demonstrate this by specific analytical calculation. 

Before continuing, we note again that the mass of the 
model belt of Ref. [jj appears to be somewhat high, 
as has been determined elsewhere 0,00- Further, it is 
known that the amount of dust is much smaller than this, 
and the gravitational mass of the Kuiper Belt is domi- 
nated by large rocks and ices. The interplanetary dust 
is actually supplied by collisions between the rocks and 
ices and lives for only of order 100,000 years in the inner 
solar system, an order of magnitude longer in the outer 
solar system. Further, the dominant mass of the rocks 
and ices is overwhelming subject to gravity and not other 
forces. Hence, there tend to be resonant concentrations 
in it vs. a smooth distribution _ E3- 

In Section ^ we will describe the gravity of spheri- 
cal balls and shells. This is followed by an introduc- 
tion to the gravity of cylindrically symmetric disks and 
rings in Section lnTl (These objects are examined in both 
the complete 3-dimensional framework and also in the 
"thin-ring" approximation, where the distribution in the 
z direction is a ^-function.) In Section llVI we apply the 
"thin- ring" approximation to both the 1 /p model and the 
Boss-Peale model pj . We then go on to full 3-dimensional 
calculations. In Sections IVl IVII and IVIII we discuss, re- 
spectively, the 1/p-density cylindrical ring, the constant- 
density cylindrical ring, and the l/r 2 -density wedge (as 
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well as the l/p 2 -density "thin ring"). We end with a dis- 
cussion where we compare the results. In particular, we 
compare the accelerations produced by the 3-dimensional 
1/p-, 1/r 2 -, and constant-density rings, as well as those 
from the Boss-Peale and 1/p 2 "thin-rings." 

We find, as expected, that neither the magnitude nor 
the shape of the Pioneer anomaly can be reproduced. 
(For comparison, in our numerical plots we will adhere 
as much as possible to the model parameters of Eq. J3Jl. 
However, since the basic formulae are analytic, they can 
be renormalized at will.) 



II. 



SPHERICAL BALLS AND SHELLS 



The 1 /r 2 gravitational force law yields that any spher- 
ically symmetric distribution with total mass M. exerts 
a force outside that distribution that is proportional to 
the total mass divided by the square of the distance to 
the center of symmetry: —GM.fr 1 . Contrarily, if the ob- 
servation point is inside a spherical distribution of mass, 
no force is exerted. 

This is an important result for understanding the ef- 
fects of a general spherically symmetric density distribu- 
tion, p(r). Since we are heading towards the 1/r distri- 
bution, consider density distributions that go as 



p{r) 



Pn(r) L n 



oo < n < oo. 



(7) 



Here {p n ,L} give the overall normalizations in terms of 
some density and length scale. These types of densities 
have long been studied by geophysicists. They often like 
to think in terms of spherical distributions and shells of 
the Earth having different functional dependences and 
thus causing different gravity signals [l2>ll3- But note: 
We are talking about spherical shells, not cylindrical 
rings. 

In the present study, we will concentrate on the dis- 
tributions for n — {0,1,2}, the constant, 1/r, and 1/r 2 
distributions. Specifically, start with the 1/r distribu- 
tion, 



pi(r) 



pi L 



It has a total mass out to a radius R of 



Mball(R) = 2n Pl L R 2 



(8) 



(9) 



(Of course, if the density distribution went to infinity 
there would be infinite mass.) From the spherical sym- 
metry condition mentioned before, we have that interior 
and exterior to the sphere 

/ . m ~G Mball{r) r< , . 

abaii{r < R) = 2 = -G2np 1 L, (10) 



abaii(r > R) 



-GMball(R) 



-G 2tt Pi L R 2 



That is, there is a constant acceleration inside the ball 
and the ordinary Newtonian inverse-square force outside 
the ball. Even so, there remains a singularity at the 
origin since there the acceleration is a non-zero constant 
pointing radially in from all directions. 

If wc now use the parameters of Ref. given in Eq. 
© above, p x = 1.74 x 10~ 16 g/cm 3 and L = 20 AU, then 
even the spherical ball of Eq. (|10fl would only produce 
an acceleration of magnitude 



abaii(r < R) 



-a 



ball 



{2irG Pl L) 



-2.18 x 10~ 8 cm/s 



(12) 



But this is smaller than ap\ So, if an entire ball of this 
density can not cause the Pioneer anomaly, how can a 
disk, let alone a ring? 

To continue, what if this were only a spherical shell 
(from i?i = 20 AU to R 2 = 100 AU)? Then, even inside 
the shell the acceleration would not be constant. By 
subtracting out the gravitational attraction of the mass 
interior to radius R\ the acceleration is 



a s heii(0 <r < Ri) 
a s heii(Ri < r < R 2 ) 



0. 



= -G 2tt Pi L 



-G 2np 1 L 



(13) 



GMball(Rl) 



G 2ir Pl L R\ 



a s heii{r > R2) 



-G 2ir Pl L (Rl - Rj) 



(14) 
(15) 



(11) 



where we write 

a s heii(r) = -(2nG piL) g shen (r) = -Cbaii g s heii(r), 

(16) 

M she ii = 2irp 1 L(R 2 -R 2 )=60M nng 

= 7.03 x 10 29 g. (17) 

Therefore, there is a constant acceleration towards the 
center of a spherical 1/r-density distribution of matter 
given by Eq. (JSJ) only if the mass distribution goes all 
the way into the origin; that is, if it is a spherical ball, 
not a spherical shell. In Figure^wc show — a s /,. e /;(r) vs r 
for the values {Ri,R 2 } = {20, 100} AU. This figure will 
be useful for comparison when we go to rings. 

Particular values of the acceleration are 

- 10 8 o,7«u({10, 60, 120} AU) = {0, 1.94, 1.45} cm/s 2 . 

(18) 

However, even here with only the first 20 AU of the 100 
AU ball deleted, the acceleration varies by an order 10% 
in the outer half of the shell and rapidly decreases to zero 
interior to that. 



3 



1.5 



0.5 




100 120 140 



We denote these various choices by: 

r y (/X) ^ z ) ' ecliptic (r, 0, 0). (23) 

(Note for future reference that, with cylindrical symme- 
try, the volume element, p, cancels the (1/p) of a pi(p) 
density function.) 

B. "Thin-ring" approximation 

1. General thin rings 

As an initial step, we start in the next section by using 
an analytic approximation, 



FIG. 1: A plot of —a s he.u{r) in units of 10 8 cm/s 2 vs. r in 
AU. 



III. CYLINDRICAL DISKS AND RINGS 

A. Full 3-d disks and rings 

Now we go on to disks and rings. We use a method 
inspired by techniques to analyze 0| cylindrically- 
symmetric objects in laboratory big-G experiments [l5|- 
|l8j. The general potential functional and acceleration 
from a cylindrical symmetric ring are 



V(r) 
V(r) 



V(r)/m test , 

/+D pR 2 
dz dp p p{p) x 
-D JRi 
2tt jj, 



(19) 



o \Jp 2 + r 2 — 2r x p cos <fr + (r z — z) 2 



a x (r) = 
a z {r) = 



-G— [V(r)] , 
dr x 

-G±Mr)h 



(20) 
(21) 

(22) 



In the above, by convention we take the component of 
the direction to the test mass in the plane of the ecliptic 
to be along the x axis: ri XjV \ — > r x . This is useful since 
we will concentrate on the case of axial symmetry. We 
also observe that the z-component of the acceleration in 
Eq. i|22fl. for gene ral positions out of the ecliptic, is easier 
to handle [ral20| in the "thin-ring" approximation of the 
next subsection. 



p(r) 2D6(z) p{p). 



(24) 



We can do this because z is generally small compared to 
p so the change in the overall result should be small and 
still symmetric about the z axis. 
This yields 



rR2 

V t hin(r) = -2GD / dppp(p) x 

JRi 



o \/p 2 + r% — 2r x p cos (f> + rf 



(25) 



athin(r) = 2GB 



dr. r 



dp p p(p) x 



o \/p 2 + r x — 2r x p cos <f> + r 2 z 



(26) 



2. Taking the v-derivative first 

One tack that can be taken (and will be in Sections 
II V Al and II V Bl below) is to first perform the r^-derivative 
in Eq. I^o]). 



a t hin{r) 



R.2 

AGD I dp p p(p) x 
Ri 

(r x — p cos < 



o [p 2 + r 2 — 2r x p cos (f> + r 2 .] 3 / 2 



(27) 



and then do the (/(-integral. Going to the plane of the 
ecliptic, r z — > 0, the result is 



rR2 

athm{r) = -4GD / dp p p(p) 

JRt 



K ( \l r 2_2 p P r+p 2 



~\J t 2 ~ 2p r + P 2 r(r + p)\J r 2 — 2pr + ; 



(28) 
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Rl 

= -4GD I dp p p{jp) 



K 



4pr 



r 2 -\-2pr-\-p [ 



4pr 



r 2 -\-2pr-\-p' 



(29) 



where the last two equalities are related by 8.127 of Ref. 
kind (see 8.113 and 8.114 in 0) are 



r(r + p) r(r — p) 

|2lj and the complete elliptic integrals of the first and second 



t. — 
1 + ~64~ 4 

"64~ 



(2n-l)!! 
2™ri! 

(2ti-1)!! 



2™n! 



i 2 " + 



p. 
In- l 



(30) 
(31) 



r 



This yields a physically intuitive p-integration that can 
be handled numerically 1]. We will use Eq. I|29l) in 
Sections HVEl and ITvTl below. 



IV. 



SPECIFIC THIN RINGS 



A. Analytic, thin-ring, 1/p-density model 



Returning to Section llll B II it turns out that the thin- 
ring, 1/p-density problem is analytically solvable. If one 
does the (^-integral before the r-differentiation in Eq. I|26() 
one can also do the second integral. (Again note, for 
this 1/p-density case, the density function cancels the 
p in the volume element, making the integrals simpler.) 
Proceeding, the potential functional in the plane of the 
ecliptic is 



d<p 



pR.2 p2n 

V r/ p(r) = -G Pl 2DL / dp / — 

JRi Jo v p + r A — Zrp cos <p 



The (j) integral is analytic and is (3.674.1 in 21] ) 

4, 



h{r > P) 
h(r<p) 



K 



T 

4, 



(?)■ 



K 

p \p 



(32) 

(33) 
(34) 



(Eqs. Q32JI and l|33(l demonstrate that for very large r the 
potential goes to —GM r i ng /r, as it should.) 

This means that the potentials outside, within, and 
inside of the ring are 



V T/P (i?2 < 



R 2 



dp 
<Ri r 

V T/p (Ri <r<R 2 ) = -8G Pl L D 

1,112 dp 



sG'/mI D j — K(^) , (35) 



r 



K 

p \p 



V T / P {r <Ri) = -8G Pl L D 



dp_ K I r 



Ri P 



(36) 
(37) 



Changing variables to t = p/r or r/p, respectively, and 
using the properties of the complete elliptic integral, the 
acceleration (ay/p = —dVT/ p /dr) is 



ar/ P (r) = -Ci 9r/p(r), 



(38) 



Ci = 8G Pl L = (4/tt) C ba ii = 2.779 x 10~ 8 cm/s 2 ,(39) 

/ n „ \ DR 2 (R 2 \ DRx (Rx\ 
9t/ p {R2 <r) = — — K — - — — K — ,(40) 



(41) 
(42) 



This acceleration is noi a constant for (i?i < r < R 2 ). 
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FIG. 2: A plot of — a,T/ p (r) in units of 10 8 cm/s 2 vs. r in 



AU. 



Putting the remaining distances in terms of AU, in 
Figure|21we plot — a*jy p (r) vs. r using the parameters of 
Ref. [5|- One can note the general features. Most impor- 
tantly, the size of the acceleration within this model of 
the Kuiper Belt is about a factor of 100 smaller than the 
anomaly! In particular, specific values of the acceleration 
are 



-a T/p ({10,60,120}AU) 
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{-0.0309, +0.0610, +0.0338} x 10~ 8 cm/s, (43) 

which can be compared to the values from a shell given 
in Eq. 1)18(1. The acceleration within the ring is of order 
40 times smaller than that within the shell. 

Observe that aT/ p (r) manifestly has other appro- 
priate physical properties. First, a-r/ p (i?2 r) — > 
—GAiring/r 2 . Next, as it should on physical grounds, 
-a T/p (r -> 0) — > 0_. Analytically, Eqs. (55)l. and 

(|42[) show that — a,T/ p {r) is slightly negative as r — > and 
goes to zero in the limit. 

One also sees the breakdowns at r = {i?2,i?i} where 
the F are singular because the arguments are unity. 
(Here and later we will cut off the heights of the 2-d 
spikes.) As we will see, these singularities result from 
having only a 2-d approximation for the non-smooth 
(hard-edged) ring. When the density is continuous in 
the p variable the spike singularity in the acceleration 
disappears, even for 2-d problems. When the problem is 
3-d, the spikes become finite cusps. (See Section jV]) 

As observed, far out a T / p (r) goes as 1/r 2 . As one 
comes in, approaches, and then passes r = i?2, the quan- 
tity —a T / p (r) starts to decrease since less mass is interior 
to the test point. Within the interior of the ring, for a 
short distance a T / p (r) is "roughly," but not exactly, flat. 
(It will be less constant in the true 3-d calculation.) Fur- 
ther, as one gets closer to R\ the acceleration changes 
sign because more mass begins pulling out rather than 
in. As predicted one sees that — a^/ p (r) is slightly nega- 
tive as r — > and it goes to zero at the origin. 



B. Another thin-ring, 1/p-density calculation 

We demonstrate here that an equivalent result for 
the 1/p-density can be obtained by the method of Sec- 
tion [TT^^I This demonstration illuminates this method 
which will be useful in the following subsection. 

If the 1/p-density given in Eq. J2J is placed in Eq. 
(|29|l . this yields the acceleration (again D will be 1 AU) 



a B p/ P {r) 



K U l ±L 



r 2 +2pr-\-p 2 



/•100 

(4GL) pi] D dp 

J2Q 



r 2 +2pr+p 2 



r(r + p) 



r(r — p) 



(44) 



= -[(8GL Pl )/2] g B p/p(r) (45) 
= -(d/2) g B p,p{r) = -(2/tt) C ball . 9sP/p (r).(46) 

The numerical integration yielding gBP/p( r ) nas to 
deal with integrable singularities at p = r, which exist 
because there the argument of K is unity. By avoiding 
the singularities, the integral is doable, except for the 
two singularities coming from the discontinuous nature 
of the ring's density at the boundaries. The result agrees 
numerically with the result in the previous subsection. 
That is, 



a-BP/pir) = a T / p (r), 



9BP/ P (r) = 2g T/p (r). (47) 



C. The Boss-Peale model 

Eq. H29|) is the integral used by Boss and Peale Q to 
study gravity from a smooth cylindrical mass distribution 
of the form 



P^ p (p-A) 2 
PBPiP) = jy 2 exp 



(p-A) 



(48) 



A = 50 AU < p < 100 AU = B, D = 1 AU. (49) 

For comparison we take this model to have the same 
mass, M r 



inng, given in Eq. (j^J. Therefore, 

-100 



Mr 



(-1UU 

AttD p% p D 2 / dp p{jp - 50) 2 

J50 

(p-AY 



exp 



4ttL> p$ p D 2 5 4 [25.8826]. 



(50) 



(If one makes the approximation that the upper limit of 
the integral goes to infinity, then the last term in the 
second line would be 26 = [r(4) + 10 T(3)].) Therefore, 



BP 

Po = 



64 



pi = (0.172) x 10~ 16 g/eni . (51) 



(25.8826) • 25 
If we place this density in Eq. I)29|l we obtain 

a.Bp(r) = -C B p gBp(r), (52) 



^BP\ 



5 3 • (25.8826) 



C B p = (4G£>p£ 

= 0.002473 Ci = (0.00687) x 10~ 8 cm/s z ,(53) 
where the quantity gsp(r) is 



gBp(r) = 



100 



dp p ip — 50) 2 exp 



(p-5Q) 



4pr 



r 2 -\-2pr-\-p 2 



E 



r(r +p) 



( ! 4 P r ~ \ 

W r 2 +2pr+p 2 J 

r(r — p) 



.(54) 



As in the last subsection, gBp{i') can be integrated 
numerically 1], but with difficulty because of the inte- 
grable singularities when r = p. The result for — asp(r) 
is shown in Figure which agrees with Figure 1 of Ref. 
|l| (except for the small, narrow spike at B = 100 - see 
below) . 

Particular values of the acceleration are 

- a BP ({10,53,73, 120}AU) = 

{-0.00686, -0.212, +0.159, +0.0325}HT 8 cm/s 2 . 

(55) 

These values, and the shape of Figure |3| reflect the dif- 
ferent type of density profile of this ring. Note that the 
curve for —a,Bp(r) is smooth when r = 50. This is be- 
cause the density varies continuously from zero at this 
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0.2 




at r = 100, which occurs since the "thin" ring abruptly 
ends there with the density Pbp{p) / P%p going discontin- 
uously from (2500 exp[— 10]) = 0.113 to zero. If the ring 
density is allowed to smoothly continue on past r = 100, 
decreasing exponentially out to infinity, the spike disap- 
pears and the resulting —a,BP 00 {' r ) becomes very slightly 
higher (lower) in magnitude than — Oflp(r) going some- 
what further out (in) from the position of the spike. 

A comparison of the normalized acceleration, app(r), 
with that for other models will be given in 

Section 



FIG. 3: A plot (solid line) of -osp(r) in units of 1CT 8 cm/s 2 
vs. r in AU. Also shown is a dashed plot of psp(p) in units 
of lCT 15 g/cm 3 . 



point. On the other hand, note the small, narrow spike 
I 



V. 3-D, CYLINDRICAL-COORDINATE, 
(l/p)-DENSITY RING 

Now we calculate the acceleration from the (1/p)- 
density in the 3-d case. Begin with the complete, exact, 
3-dimensional integral defined in Eqs. I|20(l and (|21(1 with 
the ring (l/p)-density of Eq. (J2J: 



a i/ P ( r ) = ~ G 



-d 



dr :r 



2tt 



-D 



dz 



dp p piL 



d JRi p \/p 2 +r%- 2rpcos(j)+ (z - r 2 ) 2 



(56) 



Going to the plane of the ecliptic, performing the p-integration (which is easy since the density cancels the volume 
element), and then doing the r-derivative yields 



a i/ P ( r ) = -(Ci/ 4 ) 9i/ P (r) 
ai/ P {r) ■■ 



Ci 

' 4 

Ci 
' 4 



dr 





r+D 


[ del 


) / dz x In 


. Jo 


J-D 



p — r cos 4> + \/ p 2 + r 2 + z 2 — 2pr 



COS( 



R 2 



/■7T P + D 

d(j) dz [$(r, Ra, z, 0) - *(r, Ri,z, <f>)], 

JO J-D 

cos 4>+ (R cos <j) — r)/S 



(57) 
(58) 

(59) 
(60) 

The z-integration can be done analytically using the two sets of square brackets in Eq. I|59() separately, with the 
complicated second piece adding an additional part to the first term. This yields 



$(r, R,z,, 



(R - rcoscf)) + S 



—r sin 



z 2 + r 2 sin 2 i 



Scos<p+ [-{R 2 +r 2 )coscj) + pR(l + cos 2 </>)]/5 



z 2 + r 2 sin 2 i 



S = 



y/R 2 + r 2 - 2rRcoscj) + z 2 . 



H{r, R, Z, , 



= —2 sin <b tan 



cos <p hi[Z + S] + sin <f> tan 



rZ sin ( 



S -> y / R 2 +r 2 -2rRcoscj) + Z 2 . 



R 2 + r 2 - 2Rr cos 4> + (-R + r cos </>)S 

(61) 
(62) 



Although it is technically possible to do the 0-integration, the result is so complicated that it is preferable to do 
the final integral numerically. The result, 



9i/ P (r) = / deb [H(r, R 2 ,D, 0) - H(r, R 2 , -D, 4>) - H(r, R 1 ,D, 0) + H(r, ife, —D, 
Jo 



4>)] 



(63) 



is used to obtain ai/ p (r), which is shown in Figure 01 (r — i?{i i 2}cos0) = 0. 
(The numerical singularities to be overcome occur when 
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FIG. 4: A plot of -a 1/p (r) in units of 10~ 8 cm/s 2 vs. dis- 
tance in AU, obtained from a (l/p)-density ring, with a 3- 
dimensional calculation. 



This figure again has the correct behaviour. It is a 
more delicate version of the "thin ring" shown in Fig. [21 
The most noticeable change from the "thin ring" is that 
the spikes of —ax/ p (r) near {Ri, R2] in Figure [21 become 
finite cusps at {i?i,i?2} of — a 1 / p (r) in Figure 01 The 
cusps are also less extreme compared to the spikes. This 
is because, for the 3-d ring, all the near-by mass is not at 
a point on the ecliptic, but along a line perpendicular to 
it. The proper limit can be seen by evaluating both the 



2-d and 3-d forms as r becomes large. By r = 1000 the 
two forms already agree to three significant figures. 

To summarize: A 1/p-density potential in a ring does 
not produce a constant acceleration within the ring. 

A comparison of the normalized acceleration, q 1 / p (r) 
with that for other models will be given in Section 



VI. CARTESIAN, CONSTANT-DENSITY RING 

We next consider a constant density disk. This is of in- 
terest for both physical and mathematical comparisons. 
We use cartesian coordinates because for cartesian coor- 
dinates the volume element is unity. Therefore, a con- 
stant density has the simplest integrals with these co- 
ordinates. (We already observed how the + y 2 
density cancels the \/ x 2 + y 2 volume element in cylin- 
drical coordinates.) This current calculation is similar to 
that used in Ref. [14J to study the metrology of a solid 
cylinder for big-G Cavendish experiments. 

To settle on po we take the same total mass and shape 
as the 1/p ring. This means 



2L 



Po = Pi - 



Pi/3- 



i?l + i?2 

Now proceed by using Eq. H20|) . giving 



(64) 



Vcon{r) = -Gp 



dy 



dx — 



dy 



dx 



+D 



dz 



D 



y/(x- r x ) 2 +y 2 + (z- r z ) 2 



(65) 



The two integrals represent the gravitational effect of a disk of radius R2 minus the effect of a disk of radius R\ , 
thus yielding a ring. Again, in the plane of the ecliptic (r z = 0) the acceleration is obtained by taking the negative of 
the derivative of the integrand with respect to r: 



d 



1 



dr [(r — x) 2 + y 2 + z 2 ] 1 / 2 \{r — x) 2 + y 2 + z 2 } 3 / 2 



(66) 



(Note that since one is taking the derivative of the square root of a square, one must be careful that the correct 
over-all sign emerges.) 

Now do the integral with respect to z. This yields 



Thus, 

acon(r) = -Gpo 
The x integral is 

Ir 



dy 



-R 2 



111 



I, 



z[r — x) 



1 D 



[(r - x) 2 + y 2 } [z 2 + (r - x) 2 + y 2 } 1 / 2 



-D 



dx 



dy 

Ki J-y/Rf-i 



dx 



2D(r - x) 



+D + \jD 2 + (r - x) 2 + y 1 



R-, 
Ri 



- In 



[(r - x) 2 + y 2 } [D 2 + (r - x) 2 + y 2 ] 1 / 2 ' 



R-2 



-D + ^D 2 + {r-x) 2 +y 2 



Ri 



(67) 



(68) 



(69) 
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acon(r) 
F(r,y,R,D) 



= -Gpo 



In 



Ri l-Rl 

dy F(r,y,R 2> D) - / dy F{r,y,Rx,D) 

-R 2 



-D 



D 2 + R 2 



2r v /i? 2 - y 2 } 



-D 



D 2 + R 2 +r 2 - 2r^jR 2 - y 2 } 
I 



[-D + JD 2 + R 2 + r 2 + 2r v / R 2 - y 2 } 



+D + JD 2 + R 2 + r 2 + 2r v / R 2 - y 2 } 



(70) 
.(71) 



This final integral can be done analytically usin g in - 
volved transformations similar to those used in Ref. |14| . 
But the end result is very complicated. Therefore, for 
clarity, a simple 1-dimensional numerical integral will be 
used. (As a result we leave unaddressed the implications 
of the relative sizes of r vs. {Ri,R 2 }, which implica- 
tions can play in the analytic form of this final integral.) 
We change all units to AU, e.g., change the variable y 
to t = y/D and multiply the external constants by the 
same D — 1 AU. Then, 



a-Con — —Cq 



100 



dt F(r, 4,100,1) 



100 



1-20 

dt F(r,i, 20,1) 

J -20 

-Co gcon(r), 



(72) 



Co 



Gp D = Ci/480 = 0.00579 x 10~ 8 cm/s 2 



(73) 

2 
S . 

(74) 



In Figure we show — acon(r)- Again we see the cor- 
rect general behaviour. With the 3-d calculation, the 
cusps at the discontinuous boundaries of the ring are 
large, but finite and hence physical. Interesting values 
of the acceleration are 



- a Co n({10,~ 20,60,- 100, 120} AU) x 10 8 = 
{-0.0165, -0.0870, +0.03146, +0.130, +0.371} cm/s 2 



Since the total mass is the same as for the 1/p ring, the 
acceleration should tend to the same limit as r gets large, 
and it does. 

A comparison of the normalized acceleration, acon{r), 
with that for other models will be given in Section lYlIII 



VII. WEDGE 1/r 2 (THIN- RING 1/p 2 ) DENSITY 

A. Wedge configuration 

Now we consider a wedge-shaped slice with the spher- 
ical density 



P2 L 2 



(76) 




40 60 80 100 120 140 



-0 . 05 



FIG. 5: A plot of —acon(r) in units of 10 8 cm/s 2 vs. distance 
in AU for a uniform ring. 



As before, the slice goes between R\ and i?2, except in 
spherical distance from the origin. The opening wedge 
angle is 

6> = t&n- 1 {D/R 1 ) = 0.049958 radians. (77) 
Keeping the mass of the slice the same, 
M„ ng = 2n{25)p 2 L 2 {R 2 -R 1 ), 



_ (78) 
S = sin O = l/\/40T= 0.049938, (79) 



(75) 

one has, 



p 2 = D/(5 L) Pl =(3 Pl = (1.0012) Pl . (80) 

In the plane of the ecliptic the acceleration from the 
wedge is 



ai/ r 2 (r) = — G 
Ri 



dr 
P2 L 2 



7T/2 + O p 2 T! 

d6 sin# / d(f> x 

7T/2-0Q ^o 

t 2 dt 



t 2 



\Jt 2 + r 2 — 2rt cos 4> sin 9 



(81) 



Because the density-functional again cancels the volume 
element, the t-integral yields 



It = In 



t — r sin 9 cos <j> + \/t 2 + r 2 — 2tr sin (f> cos ( 



Ri 
(82) 
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Now taking the negative of the r-derivative yields 



/■tt/2+00 

a 1 / r 2(r) = —Gp 2 L 2 / J x 

Jn/2—e 

smC 



i0 d(j> U{r,R 1 ,R 2 ,9, 
Jo 



4>), (83) 



U(r, Ri,R 2 ,0,< 



(— sin^ cos 4>) St + r — ts'mO cos< 



R 2 



[t — r sin 9 cos <j> + S t ]S t 
S t = \/t 2 + r 2 - 2rt cos sin 0. 

The ^-integral is completely analytic, and yields 



(84) 
(85) 



l<t> (r, R\,R 2l 



At 






K 


Vs~ 




it 






K 


rS+ 





>-Atr sin/ 



n R2 



Si 



Atr sin# 



si 



J R 

i?2 



(86) 
. (87) 



Ri 



S± = \/t 2 + r 2 ±2rt sin( 



We thus have 

a i/r 2 (r) = -C 2 g 1/r 2(r) (89) 
C 2 = Gp 2 L 2 /D = d/{88) = 2.5031 C x (90) 

(•7r/2+0n 

3l/r 2(r) = D d6 sin8 l (t ,(r,R u R 2 , 6). (91) 

This integral can be done numerically and is used in 
a 1 / r 2(r), shown in Figure HO (The only numerical sin- 
gularity problems are if both 9 = n/2 and also r is either 
R 2 or Rl) 
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FIG. 6: A plot of — a 1 / r 2(r) in units of 10 cm/s vs. r in 
AU. 

The most interesting observation is that this result is 
very similar to that from the 1/p-density cylindrical ring 
shown in Figurc^J (This point will be shown even better 
in Section lYl! Il l The fact that the density is falling off 



faster with distance (1/r 2 vs. l/p) is compensated for 
by the increasing spherical width, which is growing as 
r sin 9q. 

A comparison of the normalized acceleration, ai/ r 2{r), 
with that for other models will be given in Section fVIIII 



B. Thin-ring configuration 

To demonstrate the correctness of the assertion that 
the growing width of the wedge with distance caused the 
wedge to behave more like a l/p ring, we now quickly 
look at the "thin-ring" l/p 2 problem. Keeping the same 
mass as before and using the formalism of Section Iff! B 21 
yields (also see Eq. (|4"4"|) 1 



PT/pi (j>) 
PT2 

CT2 

9t/ p 2 (r) 



PT2L 2 

(R 2 - Ri) = 
Lln(R 2 /R 1 ) 

-C T 2 9T/ P 2 (r), 

(R2-R1) 



AD \n(R 2 /R 



(2.485) pa 

- - (24.85) Ci, 



(92) 

(93) 
(94) 
(95) 



D I dp 

'Ri 

E 



K 



( / 4pr ~ \ 

r 2 +2pr+p 2 J 

p r(r + p) 



Apr \ 
r 2 +2pr+p 2 J 



p r(r — p) 



(96) 



In Figure we show a plot of a,T/ p 2 {r). One clearly 
sees the difference between the 1/r 2 wedge and the l/p 2 
thin ring. With its rise going inward within the ring, 
— aT/p 2 (r) displays the higher mass concentration at r = 
Ri . (Again there is the thin-ring caveat that the spikes at 
r = {Ri,R 2 } would be finite cusps in a 3-d calculation.) 



. 2 r 
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FIG. 7: A plot of — a T/ / p 2(r) in units of 10 8 cm/s 2 vs. r in 
AU. 
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A comparison of the normalized acceleration, a-r /p 2 ( r ) - 
with that for other models is also given in Section lVlllI 

VIII. DISCUSSION 

The different physical models we have investigated in 
this paper provide an intuitive understanding about what 
type of accelerations can be obtained from Kuiper Belt 
models. In particular, they can not easily yield a con- 
stant (or even an approximately constant) gravitational 
acceleration in a cylindrical system. 

As to the specific gravitational accelerations in the 
plane of the ecliptic, a(r), we found: 

• Starting out with Figure ^ one sees that even a 
spherical shell with a 1/r density only yields an 
approximately constant acceleration near the outer 
edge of the shell. 

• Continuing on to a "thin ring" with sharp edges, 
the 1 /p density produces an acceleration that is sin- 
gular at the edges of the ring and is approximately 
constant near the middle of the ring. (See Figure 
□ 

• Contrary to this, the smoother-density, "thin-ring" 
Boss-Peale model produces a smooth acceleration 
at the inner edge and shows only a slight, narrow 
spike if the density has a small discontinuous jump 
at the outer edge instead of decreasing smoothly 
to infinity. (See Figure |3|) Thus, the physical dif- 
ferences in shape between the 1/p and Boss-Peale 
models end up being instructive. 

• The full 3-dimensional, 1/p model, yields a finite 
acceleration everywhere, so the cusps at the edges 
of the ring are finite compared to the spikes of the 2- 
dimensional "thin-ring" approximation. (See Fig- 
ure |H) 

• The 3-dimensional constant-density ring produces 
softer cusps yet a more undulatory variation than 
the 1/p ring. (See Figure 0) It is intermediate in 
its effects between the 3-d, 1/r ring and the 2-d, 
Boss-Peale "thin" ring. 

• The 3-dimcnsional wedge, with a spherical fall off 
in density of 1/r 2 , produces an acceleration that is 
very similar in shape to that from the 1/p cylin- 
drical ring. (See Figure |S|) This is because the 
growing width of the wedge with distance approx- 
imately makes up for the faster fall off of density 
with distance. 

• The above assertion is demonstrated by the con- 
trasting behaviour of the l/p 2 -density, thin ring's 
a T/ P 2 i r )i compared to the wedge's au r s{r). It 
varies much more in the belt and reaches a high 
maximum near r = R\. (See Figure 0). 



(We also mentioned how to extend these results to out 
of the plane of the ecliptic by taking r z ^ and then 
studying both a{r x ) and a(r z ) = —{d/dr z )V(r).) 
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FIG. 8: Plots of [— ffl(r)], in units of 10~ 8 cm/s 2 , vs. distance 
(in AU), for: (i) the 3-d, 1/p-density ring, [— ai/ p (r)] - short- 
dashed line, (ii) the Boss-Peale 2-d "thin ring," [-asp(r)]- 
narrow line, (iii) the constant-density ring, [—acon(r)] 
medium line, (iv) the l/r 2 -density, wedge, [— a lj/r 2(r)] 
dashed line, and (v) the l/p 2 -density, "thin ring," [— a T / p 2 (r)] 
- wide line. 

The results emphasize how difficult it is to achieve a 
truly constant acceleration within a finite cylindrically- 
symmetric system (not even considering how much mass 
would be needed to mimic the Pioneer anomaly). This 
difficulty can be put in mathematical context. Consider 
just the "thin ring," which is mathematically simpler 
than the full 3-d ring. Starting with Eq. (|25l) , a constant 
acceleration between i?i and i?2 would be produced by a 
density pc(p) that satisfied 

r = Const. / dp p p c {p) I^r), (97) 

where /^(r) is given in Eqs. (|33|l and l|34|) . That is a 
complicated inverse problem. Formally it could be done 
by a decomposition into cylindrical harmonics, but that 
is not the point here. 

Finally, in Figure |S| we show a direct comparison of 
the physical accelerations of (i) the 3-d, 1/p ring, (ii) the 
2-d, Boss-Peale "thin ring," (iii) the 3-d, constant ring, 
(iv) the 3-d, 1/r 2 wedge, and (v) the "thin," 1/p 2 ring, 
all with the same total mass, 1.96 M®. (As before we 
cut off the infinite spikes at the boundaries of the thin 
rings.) When r — > oo, all the curves tend to [GA4 r ing/r 2 ], 
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as they should. This is even though the differing density 
distributions produce quite different accelerations within 
the ring. 

To within normalizations, the results in Figure |S] agree 
with the type of results published previously for Kuiper- 
Belt disks 0, 0| • Most importantly, within the ring the 
acceleration is not constant. Further, especially in the 
central portions of the rings, the accelerations are ap- 
proximately two orders of magnitude too small to explain 



the Pioneer anomaly. 
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